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ABSTRACT 


This thesis develops the mathematical expressions 
to evaluate component weight and space for all signifi- 
cant components within the containment vessel (secondary 
shielding) of a Gas Cooled Reactor (GCR) closed cycle 
gas turbine plant for ship propulsion. The equations 
are solved for a range of operating conditions using 
helium as the working fluid. Information from these 
equations was used to obtain the following: 


1. Optimum total weight of components. 
e. Curves of component weight. 
4. Curves of component dimensions. 


The resulting data support the following conclu- 
sions: 


1. Optimum weight occurs at a pressure 
ratio that is lower than the value that gives best 
efficiency. 


eC. A heavy weight penalty must be paid 
for small increases in efficiency. 


3. The regenerator is the controlling com- 
ponent in a cycle of this type. As heat transfer re- 
quirements go up (regenerator effectiveness approaches 
100 %), regenerator weight rises sharply and limits 
appreciable gains in efficiency. 


4, Rotating machinery length will have a 
greater influence on plant design than the weight of 
this component. 


The curves should give the designer a rough 
idea of machinery component sizes and weights, and 
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are intended for use in the development and design of 
an efficient, nuclear gas turbine power plant. 


Thesis Supervisor: S. ©. Powell 


Title: Associate Professor of 
Marine Engineering 
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NOMENCLATURE 


cross section or heat transfer area (ft’ ) 


absolute velocity (ft/sec) 
constant 


specific heat at constant pressure (BTU/Lb°R) 


diameter (ft) 


pitch diameter at mid-blade 
diameter 


hydraulic diameter = S(eress section area) (4) 


wetted perimeter 


equivalent outside flow diameter per tube (ft) 


Fanning friction factor 
blade taper factor 


gravitational constant 


surface film coefficient of heat transfer (BTU/ft? -HR-°F) © 


height (ft) 
conversion factor (778 ft-Lb/BTU) 


ratio of specific heats 
constant 


thermal conductivity (BTU/ft-HR-CF) 

Heneta (LG) 

blade length 

mach number 

mass (Lb) 

number of tubes 

number of fuel rods 

number of stages 

pressure (psia) 

distance between tube centers (tube pitch), 


Fe Fa 


(ft) 





by 


Prandtl number 
heat rate (BTU/HR) 
constant, function of tube arrangement 


gas constant (ft-Lb/Lb°R) 
radius (ft) 


pressure ratio (>1) 

Reynolds number 

entropy (BTU/Lb°F) 

temperature (CR = °F + 460) and °F 


heat transfer coefficient (BTU/ft? -HR-°F) 
circumferential velocity on pitch diameter 


exponent = aa 


velocity (ft/sec) 
volume (ft’) 


exponent = al 


work (BIU/Lb) 
weight (tons) 


flow rate (Lb/HR) - 


lf 
pressure loss factor = (—t) 


5 
stress parameter = —— (ft) 
P ff 


angle of relative velocity 


blade stress (psi) 


efficiency 


regenerator effectiveness 


angular velocity (rad/sec) 
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ratio of compressor isentropic exit temperature to 


its inlet temperature 


density (Lb/ft’) or (tons/ft’ ) 
viscosity (Lb/HR-ft ) 


log mean temperature difference 


subscripts 


A 


el 


pc 


accumulator 


axial 
area 


compressor 
cold 


diameter 

gas 

high pressure compressor 
high pressure turbine 
hot 

inside 

intercooler 

low pressure compressor 
low pressure turbine 
length 

mean 

number of tubes 

outside 

piping 


precooler 





RG 


reactor 
regenerator 


surface 
stage 


turbine 
reactor vessel 


water 
weight 


cross section 


state Points 


Le 


ap 


es. 


7S. 


8s. 


precooler outlet 
LPC inlet 


LPC exit 
intercooler inlet 


isentropic LPC exit 


intercooler outlet 
HeC anlet 


HPC exit 
cold side regenerator inlet 


isentropic HPC exit 


cold side regenerator outlet 
reactor inlet 


reactor outlet 
HPT inlet 


Hel exis. 
LPT inlet 


isentropic HPT exit 


hel exer 
hot side regenerator inlet 


isentropic LPT exit 


hot side regenerator outlet 
precooler inlet 
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Ll. INTRODUCTION 


The thermodynamic potential of the closed cycle gas 
turbine power plant has been well recognized for many 
years. Early attempts to effectively utilize this poten- 
tial, however, were unsuccessful because of poor compon- 
ent efficiencies and the lack of a suitable high tempera- 
ture heat source. Technical developments of high tempera- 
ture gas turbines over the past years have enabled the 
designer to develop components with the necessary thermal 
efficiency, and the advent of the nuclear reactor has 
created the type of heat source necessary to drive the 
system. 

In the design of a ship propulsion plant, certain 
parameters are not subject to the designer’s control or 
choice. The maximum temperature in the cycle (Te) may 
be restricted because of reactor fuel cladding creep 
strength limitations or the restriction imposed by the 
high pressure turbine (HPT) blading. The minimum tempera- 
ture (T,) of the working fluid is likely to be set from 
the sea water inlet temperature and reasonable size and 
weight coolers. A fixed reactor inlet temperature (T.) 
is necessary to prevent cyclic thermal stresses from 
damaging the reactor structural material and moderator. 
In a gas cooled reactor, the inlet temperature is put at 
a relatively high value to anneal out radiation damage 


in the moderator. Compressor and turbine efficiencies 
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are dependent on the state of technical development of 


these machines. In view of the above, Tes Th» Tos and 


rotating machinery efficiencies will most likely be 
fixed for any cycle because of mechanical design or 
metallurgical restraints and size of the plant. 

The designer then has at his disposal other para- 
meters which he must fix so as to give the most economic, 
efficient or the least weight and/or space propulsion 
plant. 

The closed cycle nuclear gas turbine plant con- 
Sidered in this work consists of two principal sections: 
the reactor, which is located within the primary shield, 
and the propulsion equipment which is located within 
the secondary shielding or containment vessel along 
with the reactor. The main objective of this thesis was 
to develop a weight and space optimizing procedure for 
all significant components located within the contain- 
ment vessel. The machinery components considered as 
significant with respect to the variation in their weights 
and sizes with changes in the cycle parameters are: the 
regenerator, coolers, rotating machinery, reactor, piping, 
and accumulator. The developed expressions for component 
weights and dimensions in terms of thermodynamic proper- 
ties, gas constants and flow rate enables the designer 
to compare the relative size and weight of each component 


for different working fluids. 
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The second thesis objective was to make a weight 
and size calculation of these components using the 
proposed helium-cooled, graphite moderated, Maritime 
Gas-Cooled Reactor (MGCR) power plant as a reference. 
Solution of the developed equations for a range of 
operating conditions provided data for determining 
minimum weight and space for the regenerative, inter- 
cooled gas turbine cycle. 

No gas turbine nuclear power plant has been or 
is presently being constructed for ship propulsion. 
However, with the technological advances being made 
in nuclear reactors, high temperature materials and 
rotating machinery, it is to be expected that a ship 
power plant of this type will someday be used. The 
Shortage of component weight and space data for such 
a system will present a major problem to the ship pre=- 


liminary designer who must develop an efficient nuclear 


gas turbine plant. The developed weight and space equa= 


tions along with the system study of a particular plant 


undertaken in this thesis should assist the naval archi- 


tect and marine engineer in developing the most economi-= 


cal, efficient or the least weight and/or space gas tur= 


bine nuclear power plant. 
For any ship propulsion plant, it will be neces- 
sary for the ship designer to produce machinery compon= 


ents which give the best performance for the least 
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weight, cost, space or some other criteria. The subma- 
rine designer for instance would want to optimize his 
combination of machinery components in order to achieve 
the best performance for the least weight and space. 
The merchant ship designer on the other hand would 
direct his work toward producing a plant that gives 
the best performance for the least cost. Optimization 
of a propulsion plant can be performed for different 
reasons. Knowing the mission, size and requirements 
of the ship, the designer is able to optimize the machin- 
ery on a basis that gives the best power plant commen- 
surate with the ship’s purpose. 

In an optimization procedure, the real problem 
is to reduce the number of variables so that they can 
be conveniently handled. The ideal situation is to 
have the weight or cost expression in terms of one para= 
meter or at most two independent parameters. Computa- 
tion work and time are greatly reduced with each less 
variable. Also desirable is an expression in closed 
analytic form which will greatly simplify calculations. 
However, the functions to be differentiated will usually 
not be found in this form, and the designer must then 
resort to a graphical, numerical or some other method 
of approximate differentiation. It is important that 
careful study be given all parameters in a weight or 


cost function to determine whether or not they are 
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independent. If two parameters are dependent, the method 
of evaluating the function will be different and gener- 
ally more difficult than if all variables are independent. 
The first step in eliminating the number of vari- 
ables is to assign possible values to those quantities 
that are beyond the control of the designer. Secondly 
the concept of “balance” is used to reduce the number 
of variables. The latter step can be applied to compon- 
ents which bear a specific relationship to each other. 
For example, the task of finding the surface area of 
two different components for a given power plant might 
be reduced to a single problem of selecting a balanced 
combination of these areas. This method would be pos= 
Sible if the effects of each area on cost or plant 
weight were evidenced by their influences on the same 
working fluid property. The problem of selecting a 
balanced combination of areas is done separate from the 
optimization of the plant as a whole. All other compon= 
ents of the plant must be sensitive only to this fluid 
property and not to the various combinations of these 
two areas which establish that fluid property. Once a 
relationship between balanced components is formulated, 
the units comprising the balanced group are treated as 
a Single component with a consequent elimination of 
variables. 


In this work, a number of variables were eliminated 


-5= 





from consideration of step one. The “Dalance” concept 
was used to form a relationship between the compressors 
and the high pressure turbine. This step eliminated 
the high pressure turbine and low pressure turbine 
pressure ratios as variables. Thus once a reasonable 
value for the cycle pressure loss is specified, the 
remaining variables are flow rate (w), top pressure 
level (P), overall compressor ratio (r), regenerator 
effectiveness Cn, and temperaturesin the cycle. 

Both cycle efficiency and net work per pound of fluid 
flow are independent of top pressure as shown in Sec=- 
tion II. Pressure is then of importance only in the 
way it effects weight and space of components. All 
remaining temperatures of the cycle are dependent on r. 
For a specific power output, flow rate is a function of 
the pressure ratio. With establishment of a fixed reac- 
tor inlet temperature as previously discussed, regenera= 
tor effectiveness is also a function of pressure ratio. 
Because the reactor inlet temperature is set, the cycle 
pressure ratio must always be less than the value that 
would make the turbine exhaust temperature equal to the 
reactor inlet temperature. For equal reactor inlet and 
turbine exhaust temperatures, a Us of 100 per cent is 
required. A pressure ratio giving a turbine exhaust 
temperature lower than reactor inlet temperature would 


require a we greater than 100 per cent and hence presents 
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an impossible cycle. 

In view of the above, there are two independent 
variables for a plant of this type. All weight and size 
equations were developed in terms of these two inde=- 
pendent parameters: pressure and pressure ratio. Sub- 
sequent solution of the equations for a narrow range 
of pressure ratios indicated extreme variations in 
total weight and size of components. 

The regenerator weight and size was found to 
dominate that of all other components by several orders 
of magnitude as regenerator effectiveness approaches 
100 per cent. Im order to achieve a higher efficiency 
without having to accept the large increase in weight 
contributed by the regenerator, consideration was given 
to the idea of using an additional piece of heat exchang- 
er equipment between the regenerator and reactor inlet. 
A small quantity of high temperature fluid could be re- 
moved from the reactor to further heat the bulk of the 
fluid from the regenerator. The extraction of this 
fluid from the reactor might occur after it has passed 
up the outside of the reflector in cooling the thermal 
Shield and vessel wall and before it flows along the 
surface of the fuel rod assemblies. Another possible 
extraction point for the small quantity of high tempera- 
ture fluid might be in the LPT after partial expansion 


or at the HPT exhaust. This “pled” fluid would then be 
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recirculated to the additional heat exchanger to further 
raise the temperature of the bulk of the fluid coming 
from the regenerator. <A considerable portion of the 
regenerator duty could then be taken over by the addi- 
tional heat exchanger. The resulting lower regenerator 
effectiveness would lower regenerator weight. A study 
of such a proposal is needed to indicate if the reduc- 
tion in regenerator weight, caused by a lower regenera= 
tor effectiveness, is more than the weight increase 
caused by the additional heat exchanger (with associated 
equipment) and the greater weights of other components 
(that might be caused by the slightly larger flow rate 


required in order to keep the same power output). 
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II. PROCEDURE 


A, General 

The basic procedure followed was: 

1. Develop mathematical expressions for component 
weights and sizes in terms of state point properties 
and system constants. 

2. Determine the required system constants. 

3. Solve the equations for a range of operating 
conditions. 

Steps 1 and 2 are carried out in this section 
with the detailed procedures of the development per-= 
formed in Appendices A, B and ©. Step 3 is found in 
Appendix D. 

The plant components considered, consist of the 
regenerator, precooler, intercooler, reactor, accumula- 
tor, piping, and rotating machinery comprised of the 
low pressure compressor (LPC), high pressure compressor 
(HPC), high pressure turbine (HPT) and low pressure 
turbine (LPT). The HPT drives the HPC and LPC, and the 
LPT provides shaft power. A low pressure turbine bdy= 
pass heat sink was not considered as a plant component. 
Preliminary studies Fil sieenaentee that direct bypass-= 
ing of the low pressure turbine without an additional 


*Numbersin brackets refer to references listed in 
Appendix E. 





heat exchanger is feasible. With this system thermal 
cycling conditions appear to be within the capabilities 
of the different component materials. A system requir- 
ing additional piping, but still no bypass heat sink, 
would use a split-flow scheme. A portion of the by- 
passed fluid would be piped directly to the precooler 
and the remainder would proceed through the regenerator 
as usual. The reactor inlet temperature could be kept 
more closely regulated with this latter method. Figure 
II-1 is a flow diagram of the system indicating state 
points and some of the symbols used in the development 
of the equations. Figure II-2 is the corresponding 
temperature-entropy diagram for the cycle. 

The mathematical development that follows is based 
on the perfect gas relationships and on turbulent flow 
of the working fluid in ducts and piping. In addition, 
the velocity of fluid flow in the system is considered 
to be of the order of magnitude where static pressure 
and temperature are approximately equal to stagnation 
conditions. All gas properties are evaluated at the 
mean fluid temperature in a particular component and 
are considered constant over the range of temperatures 
encountered in that component 

For minimum compressor work in a two-stage com-= 
pression with intercooling to the low temperature of 


the cycle (T)); the pressure ratio of each compressor 
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Figure IIZ-1 


Flow Diagram for helium-cooled, graphite moderated 
Nuclear Gas Turbine Power Plant 


Figure IIf-2 
Temperature-Entropy Diagram 
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can be readily shown to equal the square root of the over-= 


all pressure ratio: 


P P 1a he 
ee - gt - (st | ae 


The above condition also implies negligible pressure 
drop of the fluid through the intercooler. 
Total compressor work can be written for equal 


compressor efficiencies as: 


acer 5 
Wo a (= a) (2 
Cc 


Total turbine work for equal efficiencies is: 
W, = 1,¢ 7 (1 - 2) (2-3) 
1 Weems Tee 


Y being the pressure loss factor in the cycle defined as: 


= Togs te. (8) a(t Ps ha)". (22 2 BY co) 
T Pa P) Ps Py 2 

t is the ratio of compressor isentropic exit temperature 

to its inlet temperature and T, is the ratio of turbine 

iniet temperature to its isentropic exit temperature. 

The factors in the last term are the fluid pressure drops 

contributed respectively by the reactor, cold side of 


regenerator, hot side of regenerator, precooler and inter- 


cooler. The pressure drop being a function of working 


= eu 





fluid, pressure, and component arrangement and geometry 
is usually estimated by the designer for each component. 
The total cycle pressure drop and hence Y, for the cycle 
or a component, can be determined. 

The net work for the cycle is simply Wy = Woe 
Using (2-2) and (2-43) the net work (Wot? per pound 


of fluid flow can be written as: 





C T et, gt (1 - $2) - & Or? = 1) (2=5) 


ase 
Qin eS, Ry. To? age Te (2-6) 
C peep senile Gaamemreapatndoas 
pl an T, qT, 


The thermal efficiency of the cycle is then: 


n = “net (2=7) 





Combining (2-5), (2-6) and (2-7): 





AM 
6 1 4 5D 
ge - ee) - & (re? - 1) 
= at Ne _ (2-8) 
Te ~ is 
ual 


For reasons discussed in Section I, Tes Ts qT, 


Ny and N, are not generally subject to the designer’s 


Silas 








control. With a reasonable value assigned to Y, the net 
work and efficiency become functions of cycle pressure 
ratio only. By setting the derivative with respect to 

t equal to zero in (2=5), the optimum value of cycle 
pressure ratio (r) is found which maximizes cycle net 
work. Similarly the optimum value of r can be found 


from (2-8) which gives maximum cycle efficiency. 


B. Heat Exchanger Equipment 

All heat exchanger equipment considered is assumed 
to behave as if of the shell and tube counterflow con= 
figuration. In the regenerator the high pressure, cold 
gas flows inside the tube, with the hotter gas flowing 
longitudinally around the tube interspace. The increased 
strength of the tubes with internal pressure would neces= 
sitate this arrangement. Water flows through the tubes 
and gas around the tubes in both coolers. This arrange= 
ment would facilitate the eemee eed periodic cleaning of ~- 
the water side of the coolers. The precooler and inter= 
cooler require large cooling water flows with low water 
velocities and low pump heads. Design factors for the 
coolers follow marine steam condenser practice. 

1. Regenerator 

For a given gas, tube size and tube arrangement, 


ye 
equations (A-30) and (A-31) can be written respectively 


# Equations carrying letter designations are from the 
corresponding appendix. 
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n => 2 a a, 
Nai = Kyg “(2q) ‘ Ng aie “ae (2=9) 
pe 
1.4 
Na.L = K Ww Ee a (teers 5 Gia yt (2-10) 
a NdL 1-1 pe Cc fal 01 
4.8 
1 [44 
where Ca ced Ale and the K’s are constants 
q O 


of proportionality. Since qT. and Th were defined as 


de ee eect 
qT = nye (211) 
Jb dk 
and =f, = -S- 44 (212) 


they may be put in terms of pressure ratio and known 


cycle temperatures. From the definition of turbine work: 
ie = Te) = Cony (Te = Te.) 
7 lies 142) 
To = Tfl-0,0 - 97) (2=13) 


Using (2-13) and (A-4), (2-11) and (2-12) can be written: 


1 


T.(2-n,) = 0 
ir = ne 2s 


iS 





1 
n,Cl- #2) 1(2-3 0) + Ten, 
T, = ae nD —-X-___-2_X (2-15) 


Substituting (2-14) and (2-15) into (2-9) and 





(2-10), cross sectional and heat transfer area variations 
become functions of presure ratio (r), nN,» W, P and 
known quantities: For a specified maximum and mini- 
mum temperature of the working fluid (T. and Ths a 
fixed reactor inlet temperature (T.) and an appropriate 
value for turbine efficiency, the area variations are 
dependent on r, FP, and w. ne becomes known once the 
pressure ratio is specified, since reactor inlet tem- 
perature is fixed. 

Regenerator weight may be written as a function 


of the heat transfer area; 


W = K, (Nd, L) + K (2-16) . 


RG 2 


where K, and K, are constants and (Nd, L) is given by 
(2-10) with (2-14) and (2-15) substituted for T, and 


The 


Aerccolers 
For a given gas, tube size and tube arrangement, 


equations (A-32) and (A-33) for the precooler become: 


: Te |? 0” z (2217) 


ES ae 





ae 
T nf Os 
Nd.L = Kyayz, w AT (Tp0) E 


where A T_ = Tg - T 


pc 


Using equations (2-13) and 





g uy 








(2-18) 


(2-19) 


(2-20) 


(A=4), Tg can be expressed as: 


7 enn 
cae Ca 
Te (1-en,-n, ten ny idee Mr ) + Nyts 
Ty = (2-21) 
1 - Ube 
AT the LMTD may be written as: 
PAN + AT, 

AS I a(AT,) = - t= IS alten i (2222) 
where AT, ig the arithmetic mean temperature difference 
and “a” is dependent on the ratio sa ae > see Figure I1-3. 

Alt 
e 
} 
| 
LAT, 
me! Moe 
Bra 4° 
Tars 


me 
—o 


4 - distance along cooler 
Pigure Jl-4 
Cooler Temperature-Length Diagram 


Ze. 





For the precooler and intercooler, a mean water tempera= 
ture Cue of 80°F is assumed, based on 90°F sea water 
inlet temperature and a 20°F water temperature rise 

vo 


through each cooler. From (2])values of ‘a range from 


1.00 to 0.710 as At; goes from 1 to 10. 
Pas 


With T., and T., specified, “a”can be approximated by: 


a = 1,227 = .0006 Ty (2-23) 


for the range of gas temperatures likely to be encoun- 
tered at the precooler inlet. 
With the above simplifications, (2-19), (2-20) and 


(2-22) may be expressed in terms of pressure ratio, ee 


and known temperatures and efficiencies by using (2-21) 
and (2+23). The resulting relations for (2-19), (2-20) 
and (2-22) may then be substituted into (2-17) and (2-18) 
giving rather long equations for the cross sectional and 
heat transfer areas in terms of flow rate, pressure, 


iemessure ratio, Nee and known quantities. 


For a given working fluid, tube size and tube ar-= 


rangement, (A=44) and (A-45) for the intercooler become: 


Nao = K et ii mp. \*? rt? O-2m) 
o 6 Na JS ( id P 
qr \ tet fe aoe 
A mics ase . 
Ndou = Kya, “(AT ) (Tid 5.8 anes 





ie 





where AT, = TT, - Ts = 7, - Ty (2-26) 
J ee aed ee eo ae 
Ore o-oo 
Pic * 2 ; 5 (2-27) 


tT, represents the temperature out of the LPC and can 


be expressed from compressor work as: 


eT | (ie ae (2-28) 


No 


Combining (2-28) with (2-26) and (2-27): 


1 


At, = Ty aa (4° = ll) (2-29) 
Cc 
fee | 
Bes ja or = ars | (2-30) 
Cc 


The LMTD corresponds to that of the precooler: 


met = aT =) (2-31) 


where“a”can be approximated by: 
a = 1.814 a = OOM Ts (2-32) 


for reasonable values of gas temperature at the inter-= 


cooler inlet. Substituting (2-28) into (2-32): 


a= 1.814 = .0014 ae (1°? 1) + 1 (2-33) 
e 


oe 





. 





Combining equations (2-30), (2+31) and (2-43): 


A 


JE 
- A.a1-.oor4 2 hr S.ayay (Poke 1) 21-540) (2-34) 
e: C 


where em was taken as 80 + 460 = S4O°R, 


By placing (2-29), (2-30) and (2-34) into both (2-24) 


and (2-25), the two intercooler area variations are 


expressed in terms of flow rate, pressure, pressure 


ratio and known temperatures and compressor efficiency. 


As for the regenerator, the precooler and inter- 


cooler weights may be expressed as functions of the 


heat transfer area. For the precooler the basic weight 


equation is: 


a 1.4 rm 
a __£ : 8 ie 
We = K, «(| ee (z) + K, (2-35) 


where AT,4T and Toc can be expressed in terms of P, 


lr; 


(2= 


n,, and known quantities from equations (2-19) through 


23). 


The weight equation for the intercooler: 


1A 4 
pawl ach, r° 
W. = K.w | : ( "ic 8 ay K, (2-46) 
m 1 


where AT, OT, and rT. . can be expressed in terms of P, 


r and known values from equations (2-29), (2-30) and (2-34). 


SEOs 





C. Reactor 


The reactor weight and size equations are developed 
by considering the reactor as another piece of heat ex- 
changer apparatus. The proposed MGCR helium cooled, 
graphite moderated reactor with a heterogeneous assembly 
of fuel rods [11], [12], [14] is used as the reference 
reactor, 

For a constant heat generation flux along the 
coolant channel and radially over the reactor, the gas 
temperature rises with a constant slope versus coolant 
channel length. The coolant wall temperature will also 


rise with the same slope from (T. + AT) to (T, + a?) 


where AT is the coolant wall temperature minus the 
temperature of the coolant. The maximum wall and cool- 
ant temperature will occur at the reactor outlet with 
the outlet coolant temperature in each channel being 
equal for uniform heat generation in the core. A con- 
stant heat generation flux is substantially true for a 
relatively small reactor core using good reflector mate- 
mer [1), £3], £15) and [18]. 


The heat transfer in the reactor can be represented 


Dive: 
Q=HA, AT (2=37) 
= = i o= A 
Q=w C(% 5) (2-378) 
with AT = maximum wall temperature minus the outlet 


meh Ie 





temperature of the coolant. H is given by (A-13): 


an Kee Cane 8 
fo. 6 eee 


(2-38) 
ne pi +8 (Pr): 


D and N represent the outside diameter of the fuel rod 
and number of fuel rods respectively. The pressure 


drop through the reactor is given by (A=-17). The term 


ems 
= (“a5 | in (A-17) can not be neglected, but is 


equal to a constant for uniform heat generation and 
fixed inlet and outlet coolant temperatures. (A-17) 
then becomes: 


+6 Ky i 
Ln ae -C, -5 (4f 5] Ca 
where for low mach numbers, static pressures and tempera= 


tures can be interchanged with stagnation properties. 


The pressure drop factor (Yp) for the reactor is: 


P V 
Yp = (38) (2-40) 


Following the procedure used in obtaining equation 
(A-22) in Appendix A, the pressure drop term for the 


reactor can be written: 


Kal ye L : 
Yp= Cz “53> M (4f 5) (2-41) 


C, being a constant determined by the fixed inlet and 


3 
Zoe 





outlet temperatures and the maximum wall temperature of 
the reactor. Using (A-24) and (A-26) the Fanning fric- 
tion factor and mach number squared are represented 


respectively by: 





DNC *| 
f= .0438 | C, (2225 
Re TS 2 Tor fe-1 1 \ 
foe- |= i (wC ) i = (2-43) 
T Eo 9) RK K PND? 
A = He : Ts = constant. Combining all constant terms: 
mR 5 
2 2 
i. (we, ) eee 
uM? = C — (2-44) 
M KR K PND? 
Substitute (2-42) and (2-44) into (2-41): 
Eee. me 
7. (s2)’ oe ep (2-45) 
R M4 K R p? yl-8p4.8 
Combine (2-47) and (2-38) and write: 
2 a6 
Cerio tance) 
aoe ta Ee (nNDL) (2-46) 


QT we Spt: Spry: 


AT as defined with equations (2-47) is a fixed value 


and so (2-46) can be written: 


q D'S (Pr) (2=47> 
(wo) ON" L 





Multiplying equation (2-45) by (2-47) and solving the 


resulting expression for ND*: 


(i) (SE) PP a 


(2-37A) can be substituted for Q to give: 


(pp) (2 T.)°? 
ND? = ic (i iy” a ee (2246) 


The above equation represents the cross sectional area 
variation of the reactor. For fixed inlet and outlet 
temperatures and a given fuel element and coolant channel 
configuration, it is dependent on the flow rate and top 
pressure level in the cycle. 

The corresponding heat transfer area (NDL) equa- 
tion is obtained by solving (2-47) for L, dividing by 
D, making the appropriate substitution for (ND? )°* and 
multiplying this result with (2-49); 


NDL = Cyp, (a) (4) (Pr): 64D w(B) (tend 


In (2-50), L is the coolant channel length or the core 
height, H. Using the bare core approximation, the rela- 
tionship between height and radius for a critical cylinder 
of minimum volume is given by reference [4] as: 


Bee 


aT R = 1.85R (2-51) 


H = 


may ihe 


25/0), 





Dis fixed for a given reactor and N, the number of fuel 
rods (coolant channels) is a function of reactor cross 
sectional area: N = CyR* . Placing this and (2-51) into 


(2-50), reactor radius cubed becomes: 


R’ = Op ay “ (4) 2) ON/g) ” 2)", Ho g(Me-M) (2-52) 


Reactor weight can be expressed as: 


W, = Con R? H (22535 


where ais the overall reactor density determined from 
a reference or similar reactor. Using (2-51) and (2-52) 


reactor weight becomes: 


ra on BPG) er AY Go] Matter eon 


For a given coolant and fuel element configuration it is: 
S _ x a 
Wp = Kp g (Te - Ts) (2-55) 


The R and H used above are for the core and do not 
represent the actual size of the reactor. The overall 
reactor radius and height, exclusive of shielding, can 
be taken as some constant multiple of the core dimen-= 
Sions for small changes in core R and H. From the MGCR 
reference reactor, the reactor vessel radius (Ry) and 
height (Hy) can be represented as: 

Ry = 1.63 (R) (2-56) 

-~25- 





Hy = 3.31(H) (2-57) 


D. Turbo Machinery 


1. Compressors 

In formulating length and weight equations for 
the turbines and compressors, it was assumed that the 
length of each component varies directly with the nun- 
ber of stages and that weight is proportional to the 
product of annulus area at the low pressure end and 
the length. Therefore by picking reasonable values 


for: efficiency, B, “2 F B, B. and 4 ) mameve ler 


(B-19) in Appendix B can be written as: 


Nipo = K 1, (2? 21) (2-58) 


LPC 
Since Ti» the minimum temperature of the cycle is fixed, 


the LPC length is: 


Lipo = Xp Gee =) (2-59) 


The annulus area given by (B-12) becomes: 


A ae WE (2-60) 


w is the flow rate, r the pressure ratio and P the top 
pressure in the cycle. 
The weight equation, after applying the condition 


stated above to (2-59) and (2-60), is: 


--6= 





Whpc = K ip (4°? =- 1) (261 ) 


Using (B-1l1) and (B-15), the tip diameter at the LPC 


low pressure end becomes: 


ce 31) a ) a 


The four corresponding equations for the HPC are 


indicated below: 





Lupo = Kg (4°? - 1) (2-63) 
W ri? 

ALpo = K. p (2-64) 

Wopg = K, SS (1°? - 1) (aa) 
W “2 0, 

Dopo 7 Fa ($3) (2-66) 


The constants for each of the above equations are 
not necessarily applicable to both the HPC and the LPC. 
Kach machine has a particular constant associated with 
Ais 

e. Turbines 

As for the compressors, reasonable values of effi-e 
ciences, B, we ‘ PB B5 and £ can be assumed or ob- 

d 
tained from the turbine manufacturer. These factors wiil 


change very little over a wide range of operating conditions 


2272 





(pressure, pressure ratio, flow rate, etc.) and so the 
Size and weight equations can be expressed similar to 
those for the compressors. 

For the HPT, (B-32) gives the number of stages 


ass 
Nepn = K Gee 2) (2-67) 


Turbine length, annulus area, weight and tip diameter 
equations become respectively upon using the basic 


relations derived in Section III of Appendix B: 


Lapp 


Ky Ch (2263) 


a wh [ye res 
Aupp = Ky i (2-69) 


i = CG =ae 


rls 


a 
where C0 = 
Tenn, 
oD => 
t == it) fab = Ctr = 10 
een 26h ae (2-70) 


rl = C(t°?-1) 4 


_x. [¥ ey (ae Ty C(r°? - 1)] = (2-71) 
nel. dd (B) = 


me eres 


For the LPT, the equations are: 


[l-n, C(r*?-1)I[¥r(1-CLt°7-1})-1] (399) 
ii ate Ce 


Lipp = Ke 


Seem 





Eel Pace ee ny, Ql - g)] (2-73) 


a 


[1-n, (1-79) (1-0, C(x 7-1) (¥1(1-0( 1° 7-1-2) 


W.-K = (2-74) 
SU Yt{1-C(1°?-1)] 
oy, 
Dipr = Ka (F) fi ~ 1,01 = Se CaP) 


The constants associated with each equation will 
be different for the HPT and LPT unless considerable 


similarity exists between the two machines. 


mere 





E. PIPING AND ACCUMULATOR 

1. General 

Both piping and the accumulator contribute signifi- 
cantly toward the total weight of the components within 
the containment vessel. In addition, the space occupied 
by the accumulator will effect the size of the contain- 
ment vessel and hence the amount of secondary shielding. 

Weight equations for piping and the accumulator 
are developed in this section. The constants associated 
with these equations were not determined due to lack of 
Suitable data from a reference design. The equations 
however indicate the trends to be expected for different 
operating conditions. 

2. Piping 

Weight of piping is considered proportional to the 


internal surface area plus a constant. 
W = K,mDL + Ky (2-76) 


For negligible heat transfer, the pressure drop term for 


the piping is expressed as: 
y,~ ‘rw 5 (2-77) 


From (A-24), the Fanning friction factor for turbulent 


flow is: 


Duc. \ -- 
f= .044 =a (2-78) 


BAO 








Mach number squared is given by (A-26) as: 


(wC_)? 








2 
wis oe LT wk (2-79) 
ae D 
P 
Substitute (2-78) and (2-79) into (2-77) and solve 
resulting expression for D. 
(yo_)°O41 21 62 
: Dp GC ):375 (F) k=) 2-80 
K., =a Cw p z ( ) 
For a given working fluid, (2-80) becomes: 
2379 
De kk. 2 (TL) °e1 (2-81) 
Pp a> 


Since temperature in equation (2-81) is raised to only 

the 0.21 power, temperature does not effect the diameter 
(weight) significantly. As a preliminary conservative 
weight estimate, the maximum temperature in the cycle is 
used with (2-81). For greater accuracy, the piping length 
could be divided into two halves, a hot and a cold section. 


Equation (2-81) can now be written as: 


9/2 
Ww ral 
D = K, aT Gin (2-82) 
Combining (2-76) and (2-82): 
2372 
W bere 
"Von (2-83) 


Considering the length of piping within the con- 


tainment vessel as a fixed quantity for similar ships, 


(2-83) becomes: 


Ae 


a = 
= 
=> -_ — 
_— — 
as 
> 
—_ 
© 





(2-84) 


The constants in (2-84) were not determined since 
suitable information from a reference design concerning 
length and weight of piping was not available. 

3. Accumulator 

The accumulator weight is expressed as a function 
of the total weight of the working fluid in the cycle 
and the allowable stress in the wall material. Consider- 
ing a long cylindrical container, the surface area, upon 


neglecting the end plates, is: 


A = 2nRL (2-85) 


The volume of the working fluid is given by: 


V = RL (2-86) 
jf =e = (2267 


where m represents the total mass of the working fluid 
within the cycle. Weight of the accumulator is simply 
the surface area times the product of wall thickness 


and material density. 


Wy = Ba eelh Oy <2 (2-88) 


Using simple hoop stress analogy, wall thickness becomes: 


We ag (2-89) 
Ta 


yer 





where P. is the external (atmospheric) pressure and P 
is the pressure of the fluid within the accumulator. 
For rapid changes of load, the pressure of the fluid 
within the accumulator must be greater than the top 
pressure of the cycle. Equation (2-87) indicates a 
high pressure is also desirable from the standpoint 

of reducing the size of the container. To facilitate 
rapid changes in plant load [10], the accumulator is 
arranged in two parts: a high and low pressure sec- 
tion. When reduction in loads occur, gas is bled from 
the HPC discharge to a low pressure accumulator or 
receiver. As the low pressure flask fills, it is 
emptied by a piston type compressor which compresses 
the gas to a pressure considerably above the top press- 
ure in the cycle. The small piston compressor discharges 
to the HP accumulator. As loads again increase, the HP 
accumulator discharges into the high pressure line 
between the HPC and the regenerator. Equating (2-86) 
and (2-87) and solving for the radius squared (R° ): 


2 mRT m 
eee Ee (2-90) 


Combining (2-88), (2-89) and (2-90), the weight of 


accumulator becomes: 


~44= 





Since the temperature of the gas to be put in the ac- 
cumulator will be of the order of 200 - 300°F (the HPC 
discharge gas temperature), the accumulator material 
need only be a carbon steel. Therefore for a given 
gas and a moderate storage temperature, equation 


(2-91) reduces to: 





P-P 
ioe Cs a? (2-92) 
A K 
Op 
Because the pressure within the accumulator will be 
several orders of magnitude larger than atmospheric 
P-P 
pressure, a <1 and (2-92) becomes: 
Pe 
= Mm i 
W, = K ¢ (2-93) 


Maintaining a constant hoop stress for a fixed cylinder 
radius means that the pressure to wall thickness ratio 
must be a constant. With these restrictions, accumulator 


weight is written: 


Wy = Ky m (2-94) 


~Zie 





F. Evaluation of Constants 

The mathematical development of the weight and 
Space equations necessitated the introduction of a num- 
ber of unknown constants. Some constants depend upon 
materials used, some are proportionality constants and 
others are constants of design. 

To solve the equations, it is necessary to evalu- 
ate each constant. Accurate evaluation of the constants 
requires a thorough investigation of similar plants, 
with note taken of the specifications, weight state- 
ments, heat balances and other information. 

The equations developed under Section II and 
Appendices A and B are general in nature and not re- 
stricted to any particular working fluid or type of 
Ship. Only in the evaluation and application of the 
constants does the problem become one of a specific 
type or design. 

Since no shipboard nuclear closed cycle gas tur- 
bine plant exists, available data for evaluating the 
constants was limited. The author relied mainly on 
information from the proposed MGCR propulsion plant 
presented in references [11], [12], [14], [15] and, 
where applicable on information from existing open 
cycle gas turbine plants. The heat balance, specifica- 
tions, and weight estimates for the helium-cooled, 


graphite moderated MGCR plant allowed satisfactory 


—45- 





evaluation of the constants for all components except 
piping and the accumulator. 

The evaluation of these constants and sample 
calculations are given respectively in Appendices © 


and D. 
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Tit, RESULES 


A. Cycle Conditions 

Original results were obtained in the form of data 
from hand calculations of the developed equations with 
reference to the helium cooled, graphite moderated MGCR 
power plant. Calculations were made by varying the 
compressor pressure ratio (r) for two different reactor 


inlet temperatures (T,): Ty = 700°F and Ty = 791 °F, 


Table III-A lists the fixed values used for both sets 


of calculations. 


Table III-A 


Fixed Cycle Conditions 


16 
LPT (net power) output =21 ,500OHP Pressure drop=— =. 66 
Working Fluid =Helium Coolant Water Inlet a 
Temperature(T,, ) =70-F 
Top Pressure (P) =1000 psia Coolant Water Outlet - 
Temperature(T,.) =90°F 
Maximum Cycle Temperature - Turbine Efficiency 
Te =1400-F (n,) =.90 
Minimum Cycle Temperature : Compressor Efficiency 
(T_ ) = T = 100°F Ci) =,86 


i 5 C 


n= 





B. Graphs 
The data is presented in graphical form in accordance 


with Table III-B. 


Table III-B 





Graphs 
PhO Figure 
Total Weight and Efficiency vs Pressure Ratio ITI-1 
Regenerator Weight vs Pressure Ratio ITI-2 


Precooler and Intercooler Weight vs Pressure Ratio IIIi-3 


Reactor Weight vs Pressure Ratio IIi-4 
Rotating Machinery Weight vs Pressure Ratio NEMS, 
Regenerator Dimensions vs Pressure Ratio Il1-6 
Precooler Dimensions vs Pressure Ratio IiI-7 
Intercooler Dimensions vs Pressure Ratio ITI-8 
Reactor Dimensions vs Pressure Ratio LIl=9 


Rotating Machinery Length and Low Pressure 
Turbine Exhaust Diameter vs Pressure Ratio ITI~10 


Flow Rate (1b/HP-HR) vs Pressure Ratio leet 


Bee 





Figure III-l 
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Figure II[I-4 
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IV. DISCUSSION OF RESULTS 


A. General 

The results obtained in this thesis are based on 
the Maritime Gas-Cooled Reactor (MGCR) plant and are ap- 
plicable only to a plant of this type. Accuracy of re- 
sults are dependent on the accuracy with which system 
constants were evaluated. 

Hand calculations were performed to solve the 
developed equations for a range of Operating conditions. 
Results are presented graphically in Section III. 

For all components, the sharp increase in weight 
or cross sectional area indicated by the curves at low 
pressure ratios is a result of the large flow rates re- 
quired at these ratios. The pressure ratios which give 
a mathematical regenerator effectiveness of 100 per cent 
are approximately 43.1 and 43.8 respectively for the high 
and low reactor inlet temperatures. <A higher pressure 


ratio than 4.1 for the cycle with T- = 791 OF would pre- 


2 
sent an impossible plant since a regenerator effective- 
ness greater than 100 %o is required. A similar state- 
ment applies to the cycle with the lower reactor inlet 


temperature. 


B. Total Weight 


Total weight and efficiency are plotted against 


pressure ratio for two different reactor inlet 
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temperatures in Figure III-l1. Minimum weight occurs at 
a pressure ratio of about 2.5 and gives an efficiency 
equal to 40.5 % for Ts = 791 °F. For a reactor inlet 
temperature of 200°F a slightly higher minimum weight 
occurs at a pressure ratio equal to 4.15 with efficiency 
again equal to 30.5 Yo. For pressure ratios greater 
than that for optimum weight, the weight rises sharply 
to infinity due to the regenerator weight. This is be- 
cause of the large increase in regenerator heat transfer 
area required as turbine exhaust temperature approaches 
reactor inlet temperature. The MGCR proposed plant is 
designed for a pressure ratio of 2.6, a value slightly 
above that for optimum weight, but giving an efficiency 
of 32 Thos, 

Figure III-l indicates that a greater gain in ef- 
ficiency is possible with small changes in pressure ratio 
above that ratio which gives optimum weight. Although 
lowering the reactor inlet temperature gives a flatter 
optimum weight, overall efficiency is significantly re- 
duced since more heat must be put into the cycle for the 
Same reactor outlet temperature. In addition, the opti- 
mum weight for T, = 700°F is larger than that fora 
higher inlet temperature because of the increased reac- 
tor weight. This greater weight is caused by the increase 
in temperature rise associated with a low inlet temperature 


as indicated by equation (2-55). The lower inlet 
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temperature would also tend to add to the complexity of 


rPeaAcGOr COnNnSUrUCGIOnN« 


C. Regenerator 


Regenerator weight reaches an optimum at a lower 
pressure ratio than that for total weight, and rapidly 
approaches infinity as nN, goes to 100 on From Figures 
III-1l and III-2, it is noted that the regenerator is the 
controlling element in total weight at high pressure 
ratios. The designer desires a high pressure ratio in 
order to achieve more efficiency, but without the large 
increase in weight that must be accepted as indicated 
by these figures. 

A 100 Yo regenerator effectiveness can never be 
realized, and in reality the designer must be willing 
to pay heavily in terms of weight and space for a regen- 
erator effectiveness as low as 90 Fo. When one con- 
Siders the weight and efficiency curve for T. = Zola 
it is worthy to note that efficiency can be increased 
5 points by going from a pressure ratio of 2.5 to 3.0. 
However, for this increase in efficiency, total weight 
is more than doubled. Undoubtedly this sizeable gain 
in efficiency could not be realized for such a weight 
penalty. 

A possible solution in obtaining part of this in- 
crease in efficiency with a less drastic increase in 


weight might be to use an additional piece of heat 
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exchanger equipment between the regenerator and reactor 
inlet. A small quantity of high temperature working 
fluid could be taken from the reactor. This extraction 
might occur after the fluid has passed up the outside 

of the reflector in cooling the thermal shield and ves- 
sel wall and before it flows along the surface of the 
fuel rod assemblies. Another extraction point for the 
small quantity of high temperature fluid might be in 

the LPT after partial expansion or at the exhaust of 

the HPT. This extracted or “plea” fluid would then be 
recirculated to the additional heat exchanger to further 
raise the temperature of the bulk of the fluid coming 
from the regenerator. A considerable portion of the 
original regenerator duty would then be taken over by 
the additional heat exchanger (pre-heater). The result-= 
ing lower regenerator effectiveness would lower regenera- 
tor weight. The heat transmission in this proposed pre- 
heater would be very high since both sides of the addi- 
tional heat exchanger will be near the top pressure in 
the cycle. With the high rate of heat transmission pos-= 
Sible under these conditions, it may be possible to keep 
the size of the pre=heater quite small. Although this 
small quantity of “plea” fluid would again enter the 
reactor at a higher temperature than the bulk of the 
gas, the thermal inertia of the reactor and the large 


quantity of fluid should prevent any harmful cyclic 
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thermal stresses. In order to produce the same power, 

a larger flow rate of fluid is required with the addi- 
tional heat exchanger. In addition, a small pump 
(compressor) would be required to raise the pressure 

of the “pled” gas to reactor inlet pressure. A study 

of such a proposal would indicate if the reduction in 
regenerator weight, caused by a lower regenerator effec- 
tiveness, is more than the weight increase caused by the 
additional heat exchanger (with associated equipment) 
and the greater weight of other components associated 
with the slightly larger flow rate. Plant complexity 
with this system may be such that would not warrant its 
use. 

Figure III-6 indicates the variation of regenera-~ 
tor length from zero to infinity as pressure ratio in- 
creases and more heat transfer area is required. The 
curves indicating length in this figure are for each 


unit of the 4 unit regenerator component. 


D. Precooler 

Figure III-4 shows for a pressure ratio of 2 and 
larger that precooler weight has a slight negative 
slope. This is to be expected since there is less and 
less heat transfer area required as the turbine exhaust 
temperature is lowered with increasing pressure ratio. 
For the lower reactor inlet temperature, the precooler 
weight increases. This effect is due to the increased 
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heat transfer area required of the precooler as regenera- 
tor effectiveness decreases. 

Precooler length rises with a decreasing pressure 
ratio as shown in Figure III-7. This would be neces- 
sary since there is little temperature drop across the 
turbines and the fluid must be essentially cooled from 
the reactor outlet temperature to the fixed low tem- 
perature of the cycle. A Significant increase in 
length is required for Ts = 700°F since the precooler 


must perform more heat transfer as ne decreases. 


E. Intercooler 

The intercooler weight, as shown by Figure III-3, 
increases Slightly at higher pressure ratios. As the 
LPC discharge temperature is raised with pressure ratio, 
more heat transfer area is required to reduce this tem- 
perature to the fixed low temperature of the cycle. 
There is no change in intercooler heat transfer area 
(weight) with variation in reactor inlet temperature. 
The intercooler must always cool the gas from the LPC 
discharge temperature to the cycle’s bottom temperature. 

Figure III-8 indicates that intercooler length in- 
creases with pressure ratio. For low ratios, little 
heat transfer area is required to bring the gas back to 
the cycle’s minimum fluid temperature. 
FF, Reactor 


Reactor weight increases sharply at low pressure 
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ratios since a bulky heat generating unit is required 
to raise the temperature of the large flow rates. 

Since reactor weight is directly proportional to flow 
rate, weight will decrease with less flow. It is also 
worthy of note that for a drop of 91°F in reacGor an— 
let temperature, a weight increase of approximately 

35 tons results at a pressure ratio of 4.0. This is 

to be expected since reactor weight as given by (2-55) 
is proportional to temperature rise to the 1.4 power. 
Reactor vessel diameter and height curves have the same 
form as the reactor weight curves. The lower inlet:tem- 
perature requires larger reactor dimensions since equa= 
tion (2-52) shows that radius cubed is proportional to 


the temperature rise in the reactor to the 1.4 power. 


G. Rotating Machinery 
As indicated by Figure III-5, rotating machinery 


has a minimum weight for a pressure ratio of slightly 
less than 2.0. At low pressure ratios, the large an- 
nulus area required causes an increase in the weight. 
At higher pressure ratios, the weight increase is not 
as sharp since length of turbo machinery has less ef- 
fect on weight than the diameter or annulus area. 
Rotating machinery length being a strong func= 
tion of the number of stages also increases sharply 
with pressure ratio. Few stages are required at low 
ratios and many stages are needed at higher pressure 
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ratios. 

Tip diameter for only the low pressure turbine 
exhaust is indicated in Figure III-10 since its dimen- 
Sions will be larger than that of the other 4 rotating 
machines. 

Although the rotating machinery weight is sig- 
nificant, the length of the unit is likely to have a 
greater influence on design of the plant. A machine 
of great length would necessitate a larger contain-=- 
ment vessel with a corresponding increase in shielding 
weight. An in-line arrangement of compressors and ture 
bines would probably result in lower total rotating 
machinery weight and fewer shaft sealing, bearing and 
support problems than an arrangement where the low pres- 
sure turbine is in a geparate casing. The latter ar- 
rangement would permit a shorter overall length by 
placing the LPT beside the HPT and compressor unit. 
This arrangement would have a slight detrimental effect 
on efficiency because of losses in piping between the 
HPT and LPT, but may be necessary if containment space 


is limited. 
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V. CONCLUSIONS 


An analysis of the results supports the following 
conclusions: 

1. Optimum weight for the reference plant occurs at a 
value of pressure ratio that is lower than the value 
that gives best efficiency. 

Efficiency at optimum weight (pressure ratio of 
25>) gis 30.5 Oo. This compares with a maximum possible 
efficiency of 45.5 % at a pressure ratio of 3.1 
(corresponding to a regenerator effectiveness of 100 Yo) 
for a reactor inlet temperature of 71°F. 

For a reactor inlet temperature of 700°F , the ef-= 
ficiency at optimum weight (pressure ratio of 3.15) is 
50.5 %o. The maximum possible efficiency is 43.23 o7o 
occurring at a pressure ratio of approximately 3.8 


(corresponding to a regenerator effectiveness of 100 o). 


2. With a plant of this type, the designer must be will- 
ing to accept a heavy weight penalty for slight effici- 
ency gains as pressure ratio is increased above the 
optimum weight value to the ratio that corresponds to a 
regenerator effectiveness of 100 %o. Table V-A indi- 
cates the changes in efficiency and total weight for 
small increments of pressure ratio above that for 


Optimum weight. 
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Table V—A 


Changes in Efficiency and Total Weight vs Pressure Ratio 


Pressure ratio increment An( Zo) A W(tons) 
Reactor inlet temperature = 791 °R 

2D — CuO erie, 5 

ee © =e. 7 Iv a 

Boal — Cad 0.8 ABS, 

2.8 - 2.9 0.8 22 

209 - 3.0 0.9 290 


Reactor inlet temperature = 700°F 


Ee> = 3.25 05 4 
5225 = 3035 o4 5 
BoD — Fed 4 IKE 
eee — Deo? 0D 16 
Bee? — 5D 04 40 
5.65 = 3.75 a, 165 


3. A lower reactor inlet temperature gives a flatter but 
migher Optimum weight. <A significant rise in reactor 


weight results from a lower reactor inlet temperature. 


4, The regenerator is the controlling component in 
this cycle. As heat transfer requirements go up (regen- 
erator effectiveness approaches 100 % ), regenerator 


weight rises sharply and limits the gain in efficiency 


«5 Om 





that is otherwise possible with a larger pressure ratio. 


5. The precooler and intercooler weights and dimensions 
experience a minimum over a wide range of possible 
operating pressure ratios and do not strongly influence 


the total weight. 


6. Length restrictions imposed by available space are 
much more likely to influence rotating machinery design 


than the weight of this component. 
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VI. RECOMMENDATIONS 


The following recommendations are made: 

1. Since most of the developed equations show that weight 
is inversely proportional to pressure (to a power), the 
tendency would be to increase pressure level to reduce 
weight and plant size. However, maximum pressure limi- 
tations would be set from the high temperature, pressure 
effects on materials. For a more complete analysis an 
investigation should be made to determine the effect of 


pressure on the reactor, piping and regenerator weight. 


2. Equations should be refined to give greater accuracy. 
Accuracy can be improved as plants of this type are de=- 
Signed and built. <A survey of these plants for compon- 
ent weight and size data will contribute to a more pre= 
cise evaluation of system constants. Such a survey 
would provide information for proper evaluation of the 


piping and accumulator constants. 


4. Shielding accounts for a major portion of total 

plant weight, with the secondary shielding weight being 
largely a function of size and arrangement of components 
within the containment vessel. A thorough study of size 
and arrangement of components is necessary to ensure a 
minimum shielding weight. Reference [9] gives a suit- 
able method for the preliminary evaluation of shielding 


weight. 
=e 





4, The regenerator weight and size dominates that of all 
other components by several orders of magnitude as regen- 
erator effectiveness approaches 100 To. An investigation 

as to the feasibility of recirculating a small quantity 

of fluid from the reactor after it has cooled the thermal 
shield or from the turbine (after partial expansion) through 
an additional heat exchanger (pre-heater), located between 
the regenerator and reactor inlet, is necessary. ‘Such a 
study would indicate if the reduction in regenerator weight, 
caused by a lower regenerator effectiveness, is more than 
the weight increase caused by the additional heat exchanger 
and the greater weight of other components associated with 


the slightly larger flow rates. 
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APPENDIX A 
HEAT EXCHANGER EQUATIONS DEVELOPMENT 


I. General 

Development of equations for the heat exchanger 
cross sectional and heat transfer area variation are 
based on reference [17]. 

Pressure drop of the gas due to entry and exit 


effect in the exchangers are considered negligible. 


Ii. Regenerator 





4 —- distance along regenerator 
Figure A-1] 
Regenerator Temperature-Length Diagram 
The rate of heat transfer from the hot fluid to 
the cold fluid for equal flow rates and heat capacities 


on each side may be written in the form: 


Q= UA,AT, (A=3) 


aes 





where AT, = Ag; = AT, is the log mean temperature dif- 


ference across the counter flow regenerator shown in Figure 


A-l. 
Regenerator effectiveness is defined as: 
T. -T T -T 
Ny = — ecm (A-4) 
8 7 Ty cS ‘< Ty 
Combining equations (A-1), (A-2), (A-3) and (A-4): 
AU=s <a = mC, ( tx ) (A-5) 
m if alee 


The overall coefficient of heat transfer is written as: 


1 én(d,/d,) , 
— ao ri Of sa al (A-6) 
Ue |) Shan Oise "7 


For a gas-to-gas heat exchanger, the thermal resis- 


tance of the tube wall may be neglected in comparison with 


the film resistances and (A-6) may be written: 


H ‘ale 
be 1 Asc | a = i (a + ge i ) (A-7) 
U" HAS” Hy n % 
A the heat transfer area inside the tube. 


where A_ = 5 
S sc 
Reference [5] expresses the heat transfer coefficient 
of a turbulent fluid as: 
Ke 8 4 
H= .023 <== (Re)'> (Pr) (A-8) 
€ 
The hydraulic diameter outside of the tubes (hot 


side) becomes 


don i dq (A+~9) 
where q for triangular tube spacing 1S: 
rp) 2 
g = 2 (Bp) 4 (A~9A) 
O 


The equivalent outside flow diameter per tube can be 


written as: 


d= dq)? Chi) 
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For the cold side, Reynolds number can be written: 


Re) , = otoce CA aeiy) 
td . Nw Cc 
an Sp 
and for the hot side: 
Sw C 


Substituting (A-11) into (A-8) and noting that 
de = ds for the cold side: 


e 
08 
"(we ye 
HO 028 “e 1.8 4 (A-13) 
Nie? Adee (Pr) ° 
i Cc 
(A-8) becomes for the hot side: 
oo 08 
Hy = ee Ky et ) (A~14) 
n°e es q(Pr),°" 


Combining equations (A-7), (A-143) and (A-14): 


Oso ot d 5S! 
si Nowa; —— : f Lal) j (A-15) 
° ° i 
028 er k 
k : 
where L = =e (Pr), <~ 1 for the temperature range 
h Gaaey 


considered. The overall heat transfer coefficient then 


becomes: . 
2O28(WC,)* Ve Ok. 


U = 2 , (A-16) 
we Sq, a 8a) Ber) eet a(- i 


The pressure drop of a gas flowing ae a tube 
with combined friction and heat transfer can be expressed 


for lowmah numbers [7] as: 
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where f is the Fanning friction factor defined as 


shear stress 

PVP /2¢, . 2 

term 2) is negligible compared to 1 and is 
2 


tT 


In the regenerator and coolers, the 


dropped. With the lowmach number likely to be encountered 
in the ducting and heat exchanger equipment, static pres- 
sures and temperatures can be substituted for stagnation 
pressures and temperatures with little loss of accuracy. 
With these simplifications, (A-17) can be integrated to 


give for the-cold side: 


2s Me L 
Lu : * =e 4 ee (A-18) 
1 
and for the hot side: 
P M, * 
ea B = —-K 4 4 f. ae (A-19) 
O 


From equation (2-4), the pressure drop factor 
(Ypq? contributed by the regenerator is: 
: fe 
Y a ( 


¥ ) 
RG Py Pe 


(A-20) 


Take antilog of (A-18) and (A-19) and combine with 
(A~20) to obtain: 


Mf M,? f 
; K | ete | Mh b) - 


Ype is a constant for the working fluid used and for 
a particular top pressure in the cycle. 
After taking the logarithm of (A-21): 
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, (Mf Mf 
i ( mee h 2 | 
Si rye ery Aer A-22 
RG is a constant since Yea as noted above, is constant. 


The relationship for the Fanning friction factor 


(f) for turbulent flow can be given as: 


(Re)° 
Substituting (A-11) and (A-12) into (A-23): 
d. Nuc 16 
Ff, = .044 | —-+——2 (A-24) 
C wa, E 


ol done, 
cy = ane (“tea )* ail 


From the mach number definition and the continuity 
relation, M* forthe cold and hot side may be written 


respectively as: 


ae 
AI 
he 2 c /{K-1 1 
C (3 Eo Cw y RK ( K star) ( ) 
fi 2 
4aj \* 4 23 Ge -1 1 
Me = ( } (wC_) é .) (A-27 ) 
h uf Eo p RK K P Nad, 
where qT. = ire . wn and Ty = Tg 7 Tg : 


Substituting (A-24), (A-25), (A-26) and (A-27) into 
(A-22) and simplifying, the following expression is 


obtained: 


£Ec= 





(wc) . 

; Neer Tt /a; 

- K-1) oe eee Gane i Wea (4) 

Yr ~ SRG (K R git 54:8 PF fo @ te ee) 


i 
where it was assumed Ce uae (Hedy at the mean 


temperatures of the cold and hot side, and all con- 
stants have been written as CRG: The pressure term 

in the denominator of (A-28) is the top pressure level 
in the cycle. 


Combining (A-3), (A-5) and (A-16);: 


yo Ux. (028 +028 (wO,)°?  k** “6 
Dp TEs @ g\tNd. yt (A-29) 
N° 8 a re eal 8p : 1 a (2 


Divide equation (A-28) by (A-29) and solve the 


resulting expression for Nd, 


“2 1.5 5c 
2 we " k= =o o) 
Nai = Cyp (7 | (5 (Pr)*" § 14a (32 
, 
fot fae 
ES G (A-30) 


Equation (A-30) represents the cross sectional area 
variation of the regenerator in terms of gas properties, 
tube size and arrangement, regenerator effectiveness, 

flow rate, top pressure level, cycle pressure ratio, and 
mean temperatures on the hot and cold side. Ong represents 


the constant terms. 
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The heat transfer area variation (Na, L) is obtained 
as a function of the same properties by solving (A-29) for 
L/d, , making the appropriate substitution for (Nai )°¢ 


and multiplying this result by (A-30): 


_ x K=1 | Rela falns 
Na,L = Cyaz wC,, | | aS } (Pr) ( *) 
? Pe lialinet 0.4 
aes ae T, /a,\t?® (A~31) 
ah ik O oe en once 
Te as 1+q a. Tite (33 
jee aie q O 


It is noted that equations (A-30) and (A-31) simplify 
considerably for a given gas, tube size and tube arrange- 
ment. The resulting expressions then become functions of 
flow rate, regenerator effectiveness, top pressure level, 
pressure ratio and mean temperatures on the hot and cold 
side. In the last two equations, the gas properties may 
be evaluated at the mean temperature on the hot or cold 
Side since Pr and the ratio of specific heats can be 
considered constant over a large range of temperatures 
for most gases. Also the change in thermal conductivity 
to the 0.2 power would vary only slightly since k under- 


goes little change with limited temperature variation. 


mil. Coolers 
Development of the corresponding cross sectional 
and heat transfer area equations was done using the 


same analogy as for the regenerator. For the coolers 
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the heat transfer coefficient was based on the outside 
of the tubes (gas side). With the approximation that 
the tube wall and water film thermal resistances are 
small compared to the gas thermal resistance, the over- 
all heat transfer coefficient was based solely on the 
gas heat transfer coefficient. 

For the precooler, the cross sectional area vari- 


ation is then given by: 
Ap \2 fp \+? “2 z 
Nd. = Cyq Ww == | (Fae fea CF2 ares (A-32) 
At R K q P 
The heat transfer area variation is: 
2 
1.4 4 les 2 eo 
Nd Lb = Gur, Ww At. “pe| fe-1) (Pr)°Ot fos on Ch=34) 
ss je R K k /(p 


Similarly for the intercooler, the following results: 


25 5 Lr y ys: 
T see ae, i, 
Na = C ie (4c) fa) iPr rt (A-34) 


O fae Are R K P 


1.4 4 1.2 2 
. 7. a : an 
= ae =] 064/10 06 7 
mo > SNaL “lar (*ss) é) ie) Pe a: ee 


where R, K, k, and Pr are the familiar gas properties. In 
the above expressions, the gas properties are evaluated at 
the mean gas temperature in the respective cooler. 


a: a TD fetes Saeed by fem ted | 
(2 bt, + 


te 





P and r are again the top pressure level and pressure 


ratio respectively of the cycle. For the precooler: 


ATM =T- 7, 


and AT, is log mean temperature difference across 
the precooler. 


For the intercooler: 
Oa TS = T, o T, = Ts, - T 


andAT is the LMTD across the intercooler. 
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APPENDIX B 
TURBO-MACHINERY EQUATIONS DEVELOPMENT 


I. General 

The assumptions involved in the development of the 
rotating machinery equations include: 

1. Uniform axial velocity over the length of the 
component. 

2. All machines of free vortex design. 

4. Equal work done per turbine or compressor stage. 

4, Static pressure and temperature equal to stag- 
nation pressure and temperature. 


5. Steady flow. 


Development of the rotating machinery equations 


are in part based on reference [6]. 


Il. Compressors 


The root stress in a blade caused by centrifugal 
force can be expressed in accordance with references 
[6] and [8] as: 

52s LEE 


= f — = 


z. (B-1) 


where f is the blade taper factor and / the density of 
blade material. Defining a stress parameter (B) as 


equal to () , (B-1) is written as: 
Si ee £ ) 
B = 26. d (4 (B-2 ) 


=) 4= 





Since the annular flow area is: 
A= nad (B~3) 
(B-2) can be written as: 
2 2 
Bee Lhe Een Us (4) (B-4) 
eng, og d 
: 18) : 
where the relation // = a72 was used to obtain the 
term on the right in (B-4). From continuity and the 
perfect gas relationship: 


WR 
a oa ‘Si 2? 


Using (B-4), (B-5) and the relation RPM =$° ae 





machinery rotational speed can be expressed: 


EG ) py) © ge 3 8° SD, 
amr 30 Pcage? Pras _— 


From (B-6) it is seen that RPM may be limited because of 
stress considerations. For a fluid with low sonic velocity, 
RPM will be limited by mach number considerations. For 
RPM designed with mach number considerations, equation 
(B-6) becomes: 
5B, 2D 
RPM = 60(~£2) a 5 (4.) (PC, )°? (B-7) 
(Ke Bt) 

after substituting (B-4) for B. All remaining turbo 
machinery relationships will be based on stress considera= 
tions. For helium as a working fluid, compressor and 
turbine design would be influenced by stress considerations 
rather than mach number limitations. 

(B-4) can be written as: 


, 


A= td (B=8) 


SA 





and d = es ” (B-9) 


Using (B-4) and (B-5), the pitch diameter is represented 


a= (2 wet qa 4) Ee g 4 ~~ (B-10) 


For the LPC inlet, byte becomes: 


iG a 
“Tae (=H) fa)” (6 4 ~ CD) 


where P and r are the top pressure level and the pressure 


Dy: 








ratio respectively of the cycle, and K represents the 
constant terms. 

Substituting (B-11) into (B-8), the inlet annulus 
area for the LPC is: 


, 7) 
i: 
hope = EF pte 877? (4) (3-12) 


For most compressors, 6, q » and PvE can be assigned 
fixed values over a range of operating conditions. Then 
pitch diameter and annulus area are dependent on the 
volumetric rate of flow and cycle pressure ratio. 

Using the same procedure, the corresponding equa-= 


marorns for the HEC are: 


2D, 

aly 25 al 2 [, b\ 7°? 
dipo = (3) (r) (ae (s a (B-13) 
= (B)7*? (£)°° (B-14) 


The tip diameter at the inlet end of the LPC can be 





al 
A = rea (rs (ay 


expressed as: r 
at eile (@), SIPC ates 


prs 





with a similar expression for the HPC. 
Applying Newton's law in the form of angular momentun, 
along with the assumptions given in section I of this ap- 


pendix, the work done per compressor stage is: 


= 
i 


U? Cc. | 
3 Bo a (tan B, - tan B,) (B-16) 


where Bo and B, are the angles as shown in Figure B-l. 


For the LPC, the work is given by: 


J OC qT 5 
Wipe = be ae aes = 1) (B-17) 


Dividing (B-17) by (B-16) the number of stages is: 


Ic T, (1°? a) 
N = Dp. nee enn nnn ee eanneanEEEE ( B-18) 
LPC Le | Ue FF 
e'se ao (tan B, - tan BL) 


where n., is the stage efficiency. (B-18) is in the 


form whereby mach number considerations might dictate 


the number of stages since 


a 
= Seo i 
ih KRg U" KRg Mv 


By substituting for U° from equation (B-4), (B-18) may 
be put in the form: 
TOMO 2) 1). (aa 
N . (B-19) 
LES au | 
1.15.8 -2 (tan By - tan B,) 
In (B-19), stress considerations now influence the number 
of stages. 
For the HPC an expression similar to (B-19) is 


applicable since the inlet temperature to both compressors 
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is equal, and only the é ratio and the angles Bo and 


RB, may vary slightly from that of the LPC. 





Compressor 





Turbine 


Figure B-1 
Stage Velocity Diagram 


By 7) 





III. Turbines 

Before developing the corresponding equations for 
the turbines, it is necessary to determine the pressure 
ratio across each turbine as a function of the overall 
cycle pressure ratio. Since the HpTf drives only the 
two compressors, the work of the HPT must equal the 


total compressor work: 
; a _ _ ae 20) 
N ppp = pC, (Te To.) = NE pee ( mt (B20) 
From equation (2-2): 


Wo = pol (7°? 1) (B-21 ) 


Equating (B-20) and (B-21) and solving for Tapp 


ee 


eel ae? ea) 


Te7EN 


where C is a constant since Ti> Tes Ny and ile will be 
fixed for a given cycle. 


From equation (2-4): 


(ty pp) (tL pp) = Ur = Yt (B=23) 
Substituting (B-22) into (B-23): 
vail = CC See (B~24) 


"LPT 

At the HP? exit (LP? inlet) the temperature can be 

expressed as a function of Te and pressure ratio. Using 
the definition of turbine work and (B-22): 


T, = 7, [1 - ae SO (B25) 


F, 


Following the procedure used for the compressors 


=7 8 





and the preceding work in this section, the pitch dia- 


meter and annulus area at the HPT exit can be written: 


WT-(1-n, (1° 2-1)] ‘f rr fy 
6 iG il 
d a B= (B-26) 
- P[1-C(1°72-1)]¥ } Cay 5 
aE 


sa 
TENEN 


wT _fi-n,C(t* 2a 1 (p)7-5(4)7 


where C = 


A = K ae 


(B-27 ) 
HPT P[1-C(1°7-1) 4 


For the low pressure turbine: 


wEC1-n, (1- pie)” I Vay pienes 
6 t Y¥ 
“LET <( E) ‘ i) ala 


w?-[l-n, (1-y4—) Ir . 9 \.5 
ee aa 1 EE VEE 
Arpp= K P : fe CB) ie (B-29) 


The tip diameter at the exhaust end of the LPT is 


given by: 2 
dD, = ( (3) ) dipm (B-30) 


For the HPT, the number of stages is: 


1 
ICpT (- ree 7 
7 6 x t ‘st 


U? Co | 
B, g (tan B. - tan B,) 


With (B-22) substituted for t (B-31) becomes: 


He 


2dCpT, C(t" 7-1) (¢} 
pac Cae B. - tan B,) 
U 





HPT = ee | 


where (B-4) was used to substitute for U’ in (B-31),C 
is again the term defined with (B-22). 


For the LPT, the number of stages become: 


-79= 











1 
so,t, (1 - 2) mere 
Nipp = ee ot (B33) 


3° | 


Using (B-4), (B-24), and (B-25): 


ZA 
ae Mp Blan oer Paprerrcnote SaNaIer— Sh a3 


IpT C 
B —F(tan Bo - tan B, )¥t[1-C(1*?-1)] 


The turbo-machinery equations developed in this and 
the preceding section appear to be quite cumbersome. How- 
ever, it should be noted that many of the factors (effi- 
Mencies 14, af, B. B. and B,) will be fixed by the tur-= 
bine and compressor design, and will be applicable to a 
wide range of operating conditions. “hen the only vari- 
ables left to the marine engineer’s discretion are Tes 


Th» t or r, Y and volumetric flow rate. 
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APPENDIX C 
EVALUATION OF CONSTANTS 


—- 


lL. General 

The MGCR reference plant has a reactor rating of 
49,7 MW (thermal) with a flow rate of 74.6 lb/sec. 
The coolant top and low temperatures in the cycle are 


respectively: 


T. = 1300°F and T, = 100°F. 
Reactor inlet temperature is Ts = 791 Re The ratio 
of turbine pressure ratio to compressor pressure ratio 
is 0.86, giving a pressure drop of 140 psi through the 
cycle. Compressor pressure ratio is 2.6 with the top 


pressure of the fluid at 1000 psia. Turbine and com- 


PeessOor efficiencies are respectively: 


de = QO.90 and A Valor 
With the above conditions, the LPT produces 20,600 HP, 
and the HPT 746 auxiliary horsepower in addition to 
driving the compressor unit. In this work, the HPT 
provides only sufficient power to overcome the com- 
pressor work, leaving the LPT with an output of 21,500 
HP for the above conditions. 

The remaining temperatures of the reference cycle 


are given below. The state points correspond to those 


indicated in Figure Il-l. 


8) 





tT, = 7, = 100°F Ta = 852°F 
T, = T, = 239°FR Ty = 301 °F 
T, = 791 °F T= 70°F 
n= 1300°R To = 90°F 
Ty = 1008°F 


Il. Regenerator 


The regenerator is composed of three equivalent 
units arranged in series. Since no heavy internal 
Supporting structures are necessary, the regenerator 
can be constructed in a U-shaped fashion to conserve 
Space when a high regenerator effectiveness is needed. 
This shape unit would also simplify the thermal expan- 
Sion problems. The constants evaluated for the regenera- 
tor dimensions and weight are based on the overall size 
and weight of the regenerator. Unit length and weight 
can then be found by dividing by 4. Each unit uses 
3/8 inch outside diameter tubes, 2OBWG, with a 5/8 
inch center=to-center triangular tube spacing. Heat 


transfer and cross section area and effective tube length 


aire < 
NdL = 23,340 ft* 
Na = 7 ft 
L = 81 ft 


Boos 





Total weight of the reference regenerator having an 
effectiveness of 0.90 is: 
Wr = 37.1 tons 
For the above tube data, q as given by equation (A-9A) 
in Appendix A is equal to 2.06. With the above infor- 
mation, equation (2-9) is written: 
>) 


Na? = K (=) i (T + .0421(r)* TP yr? Cain) 
ND “\Ty, P S ° h 


with T, and 1, defined by (2-11) and (2-12) respectively. 
Using (C-1) and the above information Kyp is evaluated 


as: 
Kyp = ieee a (C=2) 
Upon using (2-10), Kyprz, becomes: 


Kypr, = 223.5 (C4) 


Number of tubes is given by: 


cross sectional area _N@ (Coudt) 
Tw 42 T 42 
peaceeacdiets J.) m Ol. (el sp 1) 


where the denominator equals 2.45 x Or for the given 
tube data. 


Effective tube length can be found by using: 


_ heat transfer area _ NdL - 
7 tN a, - Nd, (0-5) 


Total regenerator weight as given by (2-16) is: 
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W = K, (NdL) + K 


RG 2 


where NdL is equation (2-10). Assuming headers and sup= 


port equipment weigh 4 tons per unit, K. is taken as: 


K, = 9.0 tons (C=G6) 


Using (2-16), (C=%) and the given operating con- 


ditions, K, is then evaluated as: 


al 


= Te x 107? tons/ft? (C27) 


Cross sectional and heat transfer area and the 


regenerator weight can then be written as: 


ioe 2 
wa? = .91(w)(s% ) 2 CB, + -0422(r)* B)°7? (0-8) 
x 
lies 1.4 ill > o 4 
NaL = 22435 «(5—) os (T+.0421(r) T) (C9) 
x 
Woo = 1-2) x ome cna eto (C10) 


with NadL given by (C-9). 


mit. Coolers 

Both coolers are similar in arrangement and con= 
struction. Each uses 5/8 inch outside diameter tubes, 
18 BWG, with a 7/8 inch center=to=center triangular tube 
spacing. Weight, effective tube length and cross sece= 


tional and heat transfer area for the precooler used in 


one 





the reference plant are respectively: 


seve = 14.6 tons 
ee el oeett 
Na° = 11.6 ft’ 
NaL = 7600 ft* 


Using equation (2-17), the above information and 
the conditions stated in part I of this appendix, Kyp 


is evaluated as 


Kir, = 1.48 (Ga 


ND 


With (2-18), Kypr, becomes: 


K Ths 290 (C=12) 


ND 
Number of tubes and length are given by expressions simi- 
lar to (C-4) and (C-5) respectively. 


Precooler weight is found by using (2-35): 


int \i4 ; 
is ae Dee ey Je) By 
Woe en(ae ) ( "ve (=) + Ky (2=35) 
m 
with AT, and ire given by (2-19) and (2-20) respectively. 


AT, 1Se5the leo mean temperature ditierence. 


Assuming a header weight of 5 tons: 


K, = 5 tons (Ca) 


Using (2-35), (C=-12), (C=-13) and the given conditions, 





K, is evaluated as: 


K, = 1.13 x 107? tons/ft’* 


This value of KS is of the same order of magnitude as 


that for a steam condenser [9] but somewhat smaller. 


(C-14) 


Weight, effective tube length and cross sectional 


and heat transfer area for the intercooler of the ref= 


erence design are: 


Ve = 14.4 tons 
L Se eOa ny 
de = ©isc4 Lt. 
NaL = 6690 ft* 


With (2-24) and the given conditions: 


Kyp = 3.el1 


Given conditions applied to (2-25) evaluates Kypr 23: 


Kypr = 477 


Equation (2-36) gives intercooler weight as: 


see 
ae ae 
7 ce 0 ee 
ie - Kiw oe) uke? 30° a K, 


K., for the intercooler is considered larger than that 


(C-15) 


(C-16) 


(2-36) 


for the precooler because of the higher pressure associ-= 


ated with the intercooler. 


Ko = 6.5 tons 


ae 


(C+17) 





Combining (2-36), (C-16), (C-17) and given conditions; 
Ke a@) =< 102 toneyiiy (C+18) 


Summarizing, the equations for the precooler and 


intercooler are: 


hiec OO lei. 
oy, 
Omi 
Na? = 1.48 7e | tp «2 £ Cai 
(8 ous (C-19) 
AS 
NaL = 290 wGeEp* , (C=20) 
Wo, = 1-13 x 1072 (Nan) aS (Ce) 
with NdL given by (C-20). 
Intercooler: 
oe, 
T “5 
2 g ioe De - 
Na* = 4.21 w ae | ne = (C=22) 
tr \l.4 4A 4 
NdL = 497 w ‘s| Tso ro (C=23) 
7 p-8 
m 
W, . 4 1.03 x 1072 (Nain) eens (C=24) 


with NdL given by (C=23). 


IV. Reactor 





Core radius cubed for the reference reactor was 
developed in sub-section © of Section II and can be 


written from equation (2-52) as: 





28 


R= kK, Ze (n - ae: (Gees) 
P 


With the core diameter of the MGCR reactor equal to 6.87 
ft, and for the given operating conditions, Kp is evaluated 


from (C=-25) as: 


Kp = .0176 (C=26) 


The weight of the reference ractor is approximately 


173 tons. Combining (2-51) and (2-53): 
We = fp 5-81 R’ (C-27) 


Substituting (C-25) and (C-26) into (C=-27) and 
applying the given conditions, reactor density is evalu- 


ated as: 
PR = .92 tons/(ft of reactor core) CC228)) 
Considering the reactor vessel (exclusive of 
shielding), (, = .105 tons/ft’. Substituting (C-25), 
(C=26) and (C-28) into (C=-27), reactor weight becomes: 


x 1.4 _W_ =n 
5) es 0945 (Te - tT.) pee (C=29) 


Core radius cubed is given by: 


R = ,0176 = (Pines (C=30) 


Core height can then be found using (2-51): 


H = 1.85R (2-51) 
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Reactor vessel radius (Ry) and height Hy are given 


by (2=56) randy (2 —57 J: 


Beco ae) (2-56) 


Hy = 3.31 (H) (2-57) 


V. Turbo Machinery 


In developing the turbo machinery equations, a 
large number of variables are present. Fortunately most 
of these variables are a function of the machine design 
and not a strong function of the cycle operating condi-~ 
tions. Therefore most similar machine variables will 
have values that are near those indicated below (for the 
reference plant) for a considerable range of operating 
conditions. 


Both machines: 


So A= lO 00e = la neoonmnt 


Compressor Turbine 
q (low pressure end) AS wai »20 
efficiency AA » 86 90 


Weight, machine length and tip diameter of the 
rotating machinery are considered to be of significance 
in determining the overall plant size and weight. There=- 
fore the constants associated with each of these three 


developed equations (weight, length and tip diameter) were 
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evaluated for 


each rotating machine. 


The following data apply to the machinery used in 


the reference 


inet t ) 
W (tons ) 


D (inches) 


design: 
LPC HPC ebeusl 
(oye, Owe 5 3.98 
4.8 Bic oa 5.42 
Ae ll 14.4 2042 


pine h 


2. 84 
Do DD 


Con 


The applicable length (L), weight (W) and tip dia- 


meter (D) equations from sub-section D of Section II were 


used in conjunction with the above given geometry data 


and the operating conditions to evaluate the constants. 


The following values were obtained as constants: 


For the 
L = 


LPC HPC edenih 
37 0 47.6 Tes 
Tl Ain Ul 8 156.0 
40.17 41.1 63.4 


LPC, the working equations are: 


37.0(t°7 = 1) 
mene (rt? Sa 
wr ).D 
10.7 (#) 

HPC: 


49.6(t°? = 1) 


=90= 


LPT 
Cle 


216.0 


Oreo 


(C-31) 


Ce 


(C=33) 


(C=34) 


W= 151.084 DOC) ea) 
Wee uy (3) ee 
For the HPT: 


eee Scr ues) 


w= 156 # S 
fieenoGce = 1)" 


oO 
1! 


fil = Ole 2a 
ea 


iE K 
where C = =—=—— and u = -= 
TeneNe ue 

For the LPT: 


oO 
i 


w* 721 )[1-1,0(1°7-1) 


Boat (B)°? [1-n, C(t" 7=1)] 


J 


22 
U 


27.2 fl=-n,C(x°7=1)3(¥2(1-Cft" 7-1] )-1) 


er [1-C(1° 21) ] 


216.0 
Yt[l - C(1°? ~ 1)] 


, 
wt LEAS 


Scie 


SOY 


(C36) 


(C-37) 


(C-38) 


(C-39) 


(C-40 ) 


a (o-42) 


(C-42) 





APPENDIX D 


SAMPLE CALCULATIONS 


I. General 
In order to indicate the effect of certain para- 
meters as they individually influence plant efficiency, 
weight and space, it is desirable to fix as many other 
variables as possible. As stated in Section I, the 
top and low coolant temperatures, reactor inlet tempera- 
ture, turbo-machinery efficiencies and pressure losses 
for the cycle will generally be fixed for a power plant 
of this type. From consideration of the developed equa- 
tions, the remaining parameters effecting the cycle are 
pressure (P), flow rate (w), pressure ratio (r), and 
regenerator effectiveness (n,)- However, as noted in 
Section I, regenerator effectiveness is specified once 
a pressure ratio is selected for this type cycle. 
Solution of all equations was performed for a 


LPT (net power) output of 21,500 HP with the below fixed 


conditions: 
T. * 1760°R ile GE BE 
T, Eccorr iy » 86 
T, = 1251°R 
P = 1000 psia —t = .co or Y = —20er 
tw =  70°R tyo = 90F 


BOs 





Another set of calculations were made with the 
reactor inlet temperature lowered to 1160°R, with all 
other parameters above held constant. 

Equations were evaluated for pressure ratios rang- 
ing from 1.3 to a ratio that placed the LPT exhaust 
temperature (T,) at the same level as the reactor in- 
let temperature (T.)- This condition would require a 
regenerator effectiveness (n,) of 100 Yo. Any pressure 
ratio larger than this value would present an impossible 
cycle since aie would have to be greater than 100 To. 

For equations where w and P appear, w is in 
lb/sec and P is in psia. All temperatures are in 


degrees absolute (OR). 


II. Sample Calculations 


Sample calculations are indicated for all equations 
using a pressure ratio of 2.4 and the fixed conditions 


given in Section I of this appendix. 


A. Flow rate 


Flow rate is found using equation (2-5): 


Aa 
net power = w [n, + Te, Q- 7) - “ ee &1)) 
(21,300 )550 _ ¥ 2660)(1.25) 
21.500 Sout 9001760225 0 cebyeae)~4 60)(1.251 19-11 


w = 81.1 lb/sec 
=94. 





ae Coie) (4000), ls aisneoe is 
wi 21,500 = 13-6 1b/(gp. BR) 


Boi hirc2ency 


Using (2-8), efficiency becomes 





n= 0.294 


C. Cycle Temperatures 
tT), is found using the following equation obtained 


from compressor work: 


oe 
yo tea ills: 1 Ct°~ em1)) 


T, = 560 [1 + ae Goo 1p) = eee 


rS 
i 


Tg is given by: 


= 
Tg =, (1-0, - #5 
= 1760 |1 = .90 f - cara) 
Te = 1362°R 
qT 


9 is found from a heat balance around the regnerator: 


One 





> 4 8 y 
Ty = 1362 - 1251 + 686 
Ty = 797°R 


D. Regenerator effectiveness 


n. is given by (A-4) in Appendix A: 


x 
net Seoaen 


_ 1251 = 686 
qs Tae ael eee 


Le = 0.846 


E. Regenerator weight and dimensions 


ro= 25+ 74 = 1251 + 686 = 968°R 
2 2 
dm oor 
a = -2-,—* 2 mene = 1080°R 


Putting numbers in (C-8), (C=-9) and (C-10) respectively: 


oD 
Nd* = .91(81.1) (f“B3¢ 55 (968 + .0421(2.4F 1080) 
Na’? = 5.85 ft? 


ot ee 
Nai = 22355561. 1) = (PE) oy ~8 (968+.0421(2. 4. ¥ 1080) 
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Ndi 


2 2b. 


Wee Pee eae C2720) 549 


Moe 23.8 tons 


Using (C-4) and (C-5), the number of tubes and 


length are respectively: 


N= 2:2, = 2490 tubes 
2. 55x10 
L = 12 eke, _ O ft. 





(3.14)(2490)(2) 45 


ee = = 16.67 ft/unit. 


F, Precooler weight and dimensions 


eects ||] 
SS) a ee a Veo) O 
The = 5 = 5 = 678°R 
AT, = Ty et =| 797 Ome a7. 


Ae, Pee (T) = Twi? _ 297 = 550 ~ (560 - 530) 


° Dn owe Ju 242 
A PO 
lt wl 
AT, = 103°F 


BOs 





Using equations (C-19), (C-20) and (C-21) 


aS, =5 
148) Coie) (Saé C7) 


2 


Na“ = 11.38 ft* 
eee 
NaL = 290 (81.1) ( ae (678) °4 ( 54 
Ndi = e220 ent. 
W..= 1.13 x 107? (8220) +5 
Woe = 14.3 tons 


Ths SiS) 
00461 


N = 2470 tubes 


_ 8220 
x (2470) (8) ts 


L = 20.4 dats 


G. Intercooler weight and dimensions 





Teer 
_ ee 2 686 + 560 5 
Aft, = Ts, ~ Tt) = O065= 560 == IcG 


-97= 





m iksis: 
Ln 50 
VAN ie = 20°R 


Using (C-22) through (C-24): 


)5 
Na@* = 3.21(81.1) 126)" (623) °° = 


Nd? = 14.5 ft? 


Nab = 477 (81.1) (Ey (Ee ye (2.4) 


(1000) °8 
NdL = 6590 ft° 
W., = 1.03 x 107? (6590) + 6.5 
a Sree 
N = “60461 7 2940 tubes 
ig eg 
1T( 2940 ) 8 1D 
L Pl ora ae Oe 
H. Reactor weight and dimensions 
Using (C=-29) and (C=-30): 
1.4 
Wp = ,0945 (1760 - 1251) Surya. 
(1000) °® 


~98= 





Wp = 188 tons 

Rr? = .0176 S8l-1) (1760 - 1251))*4 
(1000) ° 

R = 35,04 

R = 3,28 ft (core radius) 


From (2-51) core height is: 


Hie eros (5.26) 


H = 6.06 ft (core height) 


Reactor vessel radius and height are obtained from 


(2-56) and (2-57) respectively: 


R Lees Cages), 


y 
eee 

Ry = Ds Se iter 

Hy = 3.31 (6.06) 

Hy = 20 ft 


I. Rotating machinery 
Using equations (C=-41) through (C-42), the length, 


weight and tip diameter are given respectively for each 


ina G : 
MEG. 
L = 37 (1.19-1) 
L= 7.0 ft 
7 CBs Cent) 
We= 114 1000 (1.19-1) 


ZOO = 





4.2 tons 


| 4 
4.07 61.1 42k) 
1000 


18 inches 


47.6 (1.19 
poll ie 
Clive: 


151 000 


4.61 tons 


1000 


=) 


(2 4)2+C moe) 


0D 
41.2 ( Shad (2,4)°o? 


14.5 inches 


18.5 (1219 


Oe elites 


Sat 


158 1660 


4.2 tons 


81.1 
03.4 ( St 


= 


19-1) Lio eee NG oe 
TeseniGh OV ee 


2D 
y [1 c.822(2.1950)1 
[1 eA os 


20.7 inches 
-100= 





NS 


29 2 j1-.9(.822)(1.19-1) 10.941 (1.42) (1-.822f1.19-1]) | 
.941(1.42)(1-.822(1.19-1)] 
2.8 ft 


ee ees) 
7 1000 


b90-bereasasy IC. 9(.822)(1.19-1))[.941(1.42) (1-. 8220.19 a-.seanasa] 


941(1.42)(1-.822[1.19-1]) 


=o Ons 
5 2 
: 81.1)2.4 7 CM oe 
PSS, ( 1000 cae 8S i Saray 


D = 25.4 inches 


=~i1Q0le= 
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